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Atmospheric Effects on the Time-Varying 
Photoelectron Flux From Satellite Surfaces 


1. INTRODUCTION 


Although potentials between a satellite and the ambient medium are normally 
only a few volts (Whippte?), potentials as high as -20, 000 V during eclipse passage 
have been observed by geosynchronous satellites (Garrett), In order to analyze 
the varying potentials during eclipse entry or exit, a detailed understanding of the 
variation of the solar-generated photoelectron flux is required. As materials can 
have radically different photoemission properties as a function of wavelength, it is 
important that the amount of illumination as a function of wavelength be specified, 
In this paper, a detailed model of photoemission will be derived based on first 
principles. The model allows a determination of the photoelectron flux as a func- 
tion of wavelength. Coupled with a model of the atmospheric density, this allows 
an estimation of the total photoelectron current as a function of satellite position 
in the earth's penumbra. This latter result permits an accurate determination of 
the varying satellite potential, given the ambient particle fluxes (Ga rrett?), 

The basic theory of time-varying photoelectron flux during eclipse passage 


be i 2 Fe . 
was presented in Garrett,” Rather than repeat that development, the first section 


(Received for publication 10 June 1980) 
I. Whipple, E.C. (1965) The Equilibrium Potential of a Body in the Upper 
Atmosphere, NASA X-615-65-256,. 


2. Garrett, H.B. (1978) Effects of a Time-V ~ 
craft Potential, AFGI.-TR-78-01 


ring Photoelectron Flux on Space- 


of this paper will start with Eq. (10) of that paper, modified to include the effects 
of wavelength (an empirical model was assumed for atmospheric attenuation in the 
original study), This equation is then used to estimate the total attenuated solar 
tlux reaching ihe satellite. In the second section, the assumed cross sections for 
the atmospheric constituents will be defined, along with the atmospheric madel. 
The next two sections describe the actual computer programs utilized and compare 
their output for aluminum and tungsten with experimental results from the INJUN 5 


satellite. A listing of the computer programs is included in the Appendix. 


2. EQUATIONS AND GEOMETRICAL CONSIDERATIONS 


The geometric relationships involved in determining the attenuation of solar 
flux reaching a satellite due to atmospheric absorption and partial obscuration of 
the solar disk by the earth have been derived by Garrett. 3 The problem is simpli- 
fied if the concept of ae the minimum ray path altitude, is introduced (see Fig- 
ure 1), Xin is defined as the minimum distance from the earth that a ray of light 
attains in going fror the center of the sun to the satellite, The problem of finding 
the percentage of the so.ar disk that is obscured reduces to that of finding the 
atmospheric attenuation as a function of Xa rather than of time or position in 


orbit, which vary greatly for a given satellite. 


SAT Xm 


Xs 


Figure 1, Hlustration of the Meaning of x, the Angle Between the Sun (S) and the 
harth's (E) Center, Sy, the Distance from the Satellite (SAT) to the Center of the 
Earth, and xy the Minimum Ray Path Altitude 


Referring to Figure 2, and extending Garrett's equations to take into account 


the dependence on wavelength As it can be shown that the valu: of total flux, F 


T 


as a function of Xu is given by: 


eo Qn rs 
P(X) -f ao f ao f [F(A, X(9, 6))P(A) sin 4]da (1) 
(ay Oo AY 


where 


X(9, 0) = X(x') - R, (see Figure 2) = minimum distance of ray to arbitrary 
point on sun above the earth's surface 
x! ; gsr (cos (9) cos (x) + sin (9) sin (x) cos (0)) 


F(A, X) = solar flux as a function of A and X 


ro > angular radius of sun (function of day of year) 


P = photoemission vield for metal surface 


Hi 


wavelength 


9,0 > see Figure 2 


X = distance from satellite to center of earth. 


Figure 2, Geometric Representation of the Obscurations of the 
Solar Disk Upon Eclipse Entry. The variables are defined in 
the text 


The unattenuated photoelectron flux is given by 


wa 
F ple) = 27(1 = cos op) f F(A) P(A) da (2) 


oO 


in electrons-cm°-sec! for Fin photons-cm7*-sec7!, or upon dividing by 
-9 
6.25 X 10°? electron charges/amp-sec, F (oc) is inamp-cm ~. A useful quantity 


for comparison with observations is the percentage of photoelectron flux, given by 


F(X) 


P(X) = Fle) (3) 


The solar flux at A is obtained by multiplying the unattenuated flux (F(A)) by the 


transmission coefficient (T(A)): 


F(A) = Ee A)TA) (4) 
where 
T(A) = exp i DZ og.(A) ‘| (5) 
jel 
+a 
Ne = .(E) dé 6 
, ug n,(€) dé (6) 
-a 


number density of ‘hia species 


bes 
" 


og. = absorption cross section of jh species at wavelength » (HQ, O, © 


J O3, NO, and Ny were considered) 4% 


and € is the ray path from the sun to the satellite normal to Xin (see Figure 3). 


Thus, € = (Ry + Z) sin 9, and it is assumed that zero absorption occurs for 
Z > 400 km. 


Figure 3. Geometry for Calculation of Column 
Densities of Absorbing Species 


3. INPUT DATA 


The unattenuated solar fluxes and 0, Op», O3, and No cross sections adopted 
in this study have been obtained from compilations of various experimental data 
representative of average solar conditions by Strobel? for 1260-7500 A, Richmond! 
for 2-1027 A, and Keneshea and Huffman” for between 1027 and 1260 A. The NO 
and H,0 cross sections are from Banks and Kockarts. The solar flux values are 
shown in Figure 4, and the absorption cross sections are tabulated in the Appendix. 

For solar radiation absorption by the Og Schumann-Runge bands (1750-2050 A), 
we adopt the parameterization used by Strobel? for the transmission in each band i, 
Tr.: 


i 
Tr. = expl- f(y, 1, +6 pe (7) 
i i Oo i Oo 


3, Strobel, D.F. (1976) Parameterization of the Atmospheric Heating Rate from 
15 to 120 km Due to O9 and 03 sorption of Solar Radiation, NRL 
Memorandum Report 4 . 


4. Richmond, A.D. (1972) Numerical Model of the Equatorial Electrojet, 
AFCRL-72-0668, AD 758 196. 


5. Keneshea, T.J. and Huffman, R.E. (1972) Solar Photoionization Rate Constants 
667, AD 756480. © 


and Ultraviolet Intensities, AFCRL-72-0667, 56 ‘ 
6. Banks, P.M. and Kockarts, G. (1973) Aeronomy, Part B, Academic Press, 
N.Y. 
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and the cross section of O5 in each band: 


3 t/2 
Go a Bile (8) 
2 
where I, = column density of O54 along the ray path, and Yi 5 a, and B are 
2 


constants tabulated in subroutine ABSORB (see Appendix). 
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Figure 4. Solar Fluxes Corresponding to F = 150, and Photoemission Yield of 
Aluminum and Tungsten, as a Function of wWhYelength 


12 


Following Richmond, ¢ the solar cycle variability of the fluxes has been param- 


eterized in terms of the 10.7 cm radio flux as follows: 


P 

s 0 1 fa) n a 

Pa” Fa'Pi0.7/F io. 7, n) (9) 
where 

Fo. 7 = Fio,7 ~ 0: 0873(F 19. 7 ~ 143+) - 0.001273(F 9 7 - 143.1 
oe = flux for the nth band or line 
Fr = flux at the solar activity level represented by PG anos 
PL = a power derived from experimental data 


The photoemission yields for aluminum and dirty tungsten, also shown in Figure 4, 
are representative of experimental and theoretical data presented by Whipple, i 
The long wavelength cutoff for dirty tungsten, which is dependent on the state of 
the metal surface, and which, according to Whipple, i varies between 1°00 A and 
2600 A, is taken here to be 2300 A to best fit INJUN 5 observations (see Figure 7), 
as was empirically done by Garrett. 2 
The number densities of ©, Oo, O3, 
modified slightly within the 85 to 100 km height range to merge smoothiv with the 
U.S, Standard Atmosphere (ig7sy! above 100 km. The NO) profile is an average 


and No below 100 km are from Strobel, 


one, based on data in Swider, : 

Note that considerations of wavelengths less than 700 \ and greater than 
4500 A, as well as absorption by No, NO, and Ho, are not really germane to the 
computation of satellite photoelectron fluxes, but have nevertheless been included 
as options in the computer code for possible use in other environmental problems 


such as photoionization, atmospheric heating, and absorption by rocket effluents, 


4. FORTRAN PROGRAM 


The Appendix contains a listing of the FORTRAN program TESTL, sample 
input data, and sample output. TESTI. is a main program written for the purpose 
of illustrating the use of the working subroutines ABSORB and INTEG, Subroutine 


ABSORB computes an array of wavelength-de pendent percent atmospheric 


7. U.S. Standard Atmosphere (1976), U.S. Government Printing Office, 
fashington, D.C. 


8. Swider, W. (1972) k-region model parameters, J. Atmos, Terr, Phys. 
34:1615-1626, 
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transmission coefficients T(A) for a given minimum: ray path altitude (XS). This 
invelves interpolating the number densities of various absorbing species at altitudes 
along the ray path, finding the integrated densities of the species along the ray 
path (fq. (6)), and multiplying by the wavelength-dependent absorption cross scc~ 
tions to compute the transmission coefficient (iq. (5)). In conjunetion with 
Subroutines FUNC and ANGIN, Subroutine INTEG performs the integration i 
Eq. (1). Given a wavelength in angstroms, Subroutine PYIEI.D gives the photo- 
emission yields of aluminum and dirty tungsten by interpolating in tables. 

Program TESTU. illustrates the use of these subroutines to Compute phato- 
electron fluxes for the INJUN 5 satellite. First, the wavelengths, solar flux 


1 : . 2 
), the absorption cross sections (cm”), and number 


values (ergs-cm72-sec™ 
densities (em7?) are read in from a data file (see Appendix for sample data files), 
An option is included to replace NO with H,0 as the fifth absorbing constituent. 
(This option has been utilized to examine the possible effects of water discharges 
from repeated launches of Heavy Lift Launch Vehicles (HILV's) connected with 
proposed future Satellite Power System (SPS) activities.) The fluxes are then 
adjusted to the appropriate solar activity level, using Eq. (9), and converted from 
ergs-cm 2-sec! to photuns-cm~2-see"?, The unattenuated photoelectron current 
is computed and printed cut in amp-cm~2, Finally, by calling ABSORB and 
INTEG, P(X), the percentage of photoelectron flux as a function of Xin tS com- 
puted and printed out. Also printed are the transmission coefficients as a function 


of 2 and wavelength, and the quantity 


r 


fr (A) P(A) da 
wo 


ae 
PCURI = ——-——___—_ (10) 
be) 


ri F (A)P(A) dA 
Pa 


el 


which illustrates the relative contributions of various wavelengths to the photoelec- 


tron flux. See Appendix for a sample computer printout, 


5. SAMPLE RESULTS 


Some insight into the process of photoemission from a satellite surface can be 
obtained by cxamining results from the present model as a function of surface 
material, wavelength, and atmospheric absorbing Species. In Figure 5, the per- 


cent of unattenuated photoelectron current, given by Eq, (10), is plotted as a 
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Figure 5. Percent of Unattenuated Photoelectron Current 
for Tungsten and Aluminum, Assuming Absorption by Q, 
Oo, O3, and No, and an Fio 7 Value of 150 


function of wavelength for tungsten and aluminum, assuming Pio, 7 * 150. Here we 
see that 50 percent of the unattenuated photoelectron current is due to wavelengths 
less than 2000 A for tungsten and 3000 A for aluminum, and that 90 percent of the 
current carl be attributed to wavelengths less than 2300 A and 3600 A respectively, 
for tungsten and aluminum. The total unattenuated photoelectron currents corre- 
g amp-om"" and 0.167 x 10°78 
10.7 > 30), and 0. 305 x 107° and 
10.7 ° 250). Thus, we see that 


photoemission from a tungsten surface, due to its greater dependence on shorter 


sponding to aluminum and tungsten are 0, 293 x 107 
-2 : Pate 
amp-cm “, respectively, at solar minimum (F 


0.228 © 1078 amp-em 72 at solar maximum (F 


Wavelengths, ¢somere vattable with $ol@r aetuyigy. Since qiost of the sti e spheorie 
absorptian ceetrs below Lo@ bri altitude, the solar cvele variation of therriosphe ri 
comapostiton has been ocutted in the present study, 

he oabove differences inthe wavelength dependence of photuetisston rates for 
different iatertals suggest differences would exist wrth regard te the relative 
lrnportanee of abserloing species, boroaunstanee, from Input Table Loin the Appendix 
wee see thet €1 © aed ¢ 4 aire Strom abserbers inthe EUV, UV, and visible, 


respectively,  f ether: from the height dependences of these constituents, we 
‘sould expect that ROV, UV, and vistble radiations are abserbed at progressively 
lower heights inthe atmosphere, Phese pots are ilustrated in Figure 4, where 
tts demonstrated that Cy absorption of near-visible radiation can be neglected 
for tungsten photoemission in the case of EUV and aluminum. Note alse that the 
50 percent and 4O percent levels of unattenuated photoelectron current are attained 
above 455 km and 85 km, respectively, for aluminum, but above &5 km and 225 kin, 
respectively, for tungsten, 

In Figure 7, the percentage of photoelectron current measured by the INJUN 5 
satellite shows good agreement with computations from the present model, assum- 
ing a tungsten surface and Foo. 7 150. The only adjustment made to optimize this 


fit was to choose a long wavelength cutoff of 2300 A for the photoemission vield 


ALUMINUM TUNGSTEN 
° ie} 
—— 0,02 — — 0,02 


—-— 0,%,0, 


Xm (KM) 


° 20 40 60 80 100 O 20 40 60 80 100 
% PHOTOELECTRON CURRENT 


Figure 6, Percent of Photoelectron Current vs X),, for Aluminum and Tungsten for 
Different Combinations of Absorbing Species and an ae) 7 Value of 150 
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curve, According to an assessment of various experimental and theoretical data 
by Whipple, : this long wavelength cutoff, which is dependent on ‘he state of the 


metal surface, varies between 1900 A and 2690 A for dirty tu ¢ on. 


fe) 20 40 60 80 100 
% PHOTOELECTRON CURRENT 
Figure 7, Theoretical Calculation of 
Percent of Photoelectron Current vs 


Xm for Tungsten and an I'y9, 7 Value of 
150 Compared to INJUN 5 Observations 


6. CONCLUSION 


The model of atmospheric attenuation presented in this report, when coupled 
with models of the photoelectron emission characteristics of satellite material, 
provides an important step forward in our ability to study spacecraft charging. A 
number of increasingly more sophisticated models of the satellite charging process 
have been developed in the last year, One of those codes, the su-called NASCAP 
(NASA Spacecraft Charging Analyzer Program) model, is capable of modeling the 
complex surfaces and materials found on a typical satellite. Coupled with the 
knowledge of the solar flux as a function of satellite position, modeling of the 
penumbral passage of a satellite such as the AF/NASA P78-2 is possible-the simple 
model of Garrett? has already been used for this purpose. As the P78-2 instrumen- 


tation provides measurements of the potentials on several different types of 


17 


materials on the satellite, the model introduced here should prove extremely use- 
ful in allowing NASCAP to predict surface potentials during eclipse passage. It 
is to this end that the report has been prepared. 
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Main Program TESTI 
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cot 


rales} 


201? 


tNIA 


maha 


im} 


tay 


PRIGTAM FESTLEINPUT, OUTFIT) 
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FORMAT ESS/9 

PRINT Kote 

PRINT 4C0L 

my 292d 1-145 

FRINT QO, TACT GL 2CIV SF TCDD, CAREC DED: Sede dN 
PRINT AF AQ 

PeINT ADL 
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POINT ADDL 

ON 2982 15945135 
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PRINT AATE 

PRINT ByOt 
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PRINT ALAN 

PRINT BODE 
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PPINT AeMP 
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TF INH20=4 DEF ENE WATE? at. ANT ANSORF CDERE HEE 
IFCINH206C Tot! GO TI 106 
AAP ES,ILLART Leah $7 
FE ar EP MV PL TEA (6 


HO 103 set, 7Mar 
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TROP 6GE, ALG D ZNES SSE POVTOTENAL, 2ZBPSINE2,Z 3d FINES, 25D) 
£ CPM te SPAIN ED 20? 
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COMT TNE 


AOSUST FLEX TO SILAP ACTIVITY LEVEL, 
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AND COMPUTE UNATTIFNUATEO FHUTOFECCTRON FLUR EPCURID. 
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Main Program TESTL. (Cont.) 
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Input Table 1. Solar Flux anda as Functions of 
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Output Table 2, Absorption Coefficients as Functions of 


Wavelength 
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output Table 2.0 Absorption Coefficients as Functions of 
Wavelength (Cont. ) 
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